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The Lagrangian average (LA) of the ideal fluid equations preserves their fundamental transport 
structure. This transport structure is responsible for the Kelvin circulation theorem of the LA flow 
and, hence, for its potential vorticity convection and its helicity conservation. 

Lagrangian averaging also preserves the Euler-Poincare (EP) variational framework that implies 
the LA fluid equations. This is expressed in the Lagrangian-averaged Euler-Poincare (LAEP) theo- 
rem proven here and illustrated for the Lagrangian average Euler (LAE) equations. 
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Introduction. In turbulence, in climate modeling and 
in other multiscale fluids problems, a major challenge is 
"scale-up." This is the challenge of deriving models that 
correctly capture the mean, or large scale flow - including 
the influence on it of the rapid, or small scale dynamics. 

Meteorology and oceanography must deal with averag- 
ing in either the Eulerian, or the Lagrangian fluid specifi- 
cation. Determining the relation between averaged quan- 
tities in these two specifications is one of the classical 
problems in the physics of fluids. 

The GLM (Generalized Lagrangian Mean) equations 
of Andrews and Mclntyre systematize the approach 
to Lagrangian fluid modeling by introducing a slow + 
fast decomposition of the Lagrangian particle trajectory 
in general form. In these (exact, but not closed) equa- 
tions, the Lagrangian mean of a fluid quantity evaluated 
at the mean particle position is related to its Eulerian 
mean, evaluated at the current particle position. Thus, 
the GLM equations are expressed directly in the Eulerian 
representation. The Lagrangian mean has the advantage 
of preserving the fundamental transport structure of 
fluid dynamics. For example, the Lagrangian mean com- 
mutes with the scalar advection operator and it preserves 
the Kelvin circulation property of the fluid motion equa- 
tion. However, the Lagrangian mean also has two main 
disadvantages: it is history dependent and it does not 
commute with the spatial gradient. 

In this paper, we place Lagrangian averaged (LA) fluid 
equations such as the GLM equations into the Euler- 
Poincare (EP) framework of constrained variational prin- 
ciples. This demonstrates the variational reduction prop- 
erty of the Lagrangian mean, encapsulated in the LAEP 
Theorem proven here: 



Lagrangian Averaged Euler-Poincare Theorem 

Lagrangian averaging preserves the variational struc- 
ture of the Euler-Poincare framework for fluids. 



According to this theorem, the Lagrangian mean's 
preservation of the fundamental transport structure of 



fluid dynamics also extends to preserving its EP varia- 
tional structure ^-^]. This preservation of variational 
structure is not possible with the Eulerian mean. The 
Eulerian mean also does not preserve the transport struc- 
ture of fluid mechanics. 

The LAEP Theorem puts the approach using LA 
Hamilton's principles and LA equations onto equal foot- 
ing. This is quite a bonus for both approaches to mod- 
eling fluids. According to the LAEP Theorem, the LA 
Hamilton's principle produces dynamics that is guaran- 
teed to be verified directly by averaging the original equa- 
tions, and the LA equations inherit the conservation laws 
that are available from the symmetries of Hamilton's 
principle for fluids. 

Outline of the paper. We begin by briefly review- 
ing the GLM theory of Andrews and Mclntyre jljj. We 
then state and prove the LAEP theorem, following the 
EP framework of ||-^| . We illustrate the LAEP theorem 
by applying it to incompressible ideal fluids. Finally, we 
discuss recent progress toward closure of these equations 
as models of fluid turbulence. 

GLM (Generalized Lagrangian Mean) theory. The 

GLM theory 111 begins by assuming that the Lagrange- 
to-Euler map factorizes as a product of diffcomorphisms, 
g(t) = 3 (i) • g(t) . Moreover, the first factor g(t) arises 
from an averaging process denoted as g(t) = 3 (t) ■ g{t) — 
g(t) . The averaging process ( • ) can be reasonably arbi- 
trary, except that it must satisfy the projection prop- 
erty, so that g(t) = g(t) = g(t). Hence, a fluid parcel 
labeled by xo has current position, 

x J (x ,i) = g(t)-x =3(t) - (<K*)-xo) = 3(x(xo,t),i) , 

and it has mean position x(xo, t) = g(t)-Xo . 
GLM velocities and advective derivatives. The 
composition of maps g(t) = 5 (t) ■ g(t) yields via the 
chain rule the following velocity relation, 



(x , t) = g{t) -x = 3 (t) -x + TE ■ (§(t) -x ) . (1) 
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By invertibility, xo = g _1 (t)-x^ = <7 _1 (i)-x. Hence, one 
may define the fluid parcel velocity at the current posi- 
tion u(x^,t) in terms of a vector field evaluated at the 
mean position (x, t) as 



u{^,t)=g-g- L {t)-^ =ff.§- 1 (t).x = u t (x,t). (2) 



A - 



The velocity relation (|lj) then implies, 

u«(x,<) = ^(x,i) + ^-u L (x,i). 



(3) 



This is a standard velocity relation from GLM theory, in 
which the Lagrangian mean velocity u L is defined as 



u L (x,i) = u«(x,i) = gg- 1 (t)-x = 'g{t)g{t)- x -x. (4) 

In the third equality one invokes the projection property 
of the averaging process and finds g — g = g from equa- 
tion (0), so that u L (x, t) = ^(t)g{t)~ x = u(x,i) . Thus, 
the Lagrangian mean velocity u L coincides with u, the 
tangent to the mean motion associated with g(t). Hence, 
one may write equation (||) in terms of the LA material 
time derivative D L /Dt as 



uf(x,f) = (^+u i -V)E,x./i 



El 

Dt 



S(x,t). (5) 



For any other fluid quantity x one may similarly define 
as the composition of functions x^(x, t) = x(x^,i) = 
x(H (x, t),t). Taking the LA material time derivative of 
X^ and using the definition of D L / Dt in equation (|^) 
yields the advective derivative relation, 



(6) 



so D L x^/Dt = (Dx/Dt)t. As in equation (@) for the 
velocity, the Lagrangian mean x L of a fluid quantity 
X is defined as 



X L (x,i) = X «(x,<) - X (x«,t) = x(ff(*) ' xo,t) • (7) 

Taking the Lagrangian mean of equation (^) and again 
using its projection property yields % L = D L x L /Dt = 

(Dx/Dt) = x . Thus, the La grangian mean defined in 
(0) commutes with the material derivative. 
Transformation factors of advected quantities. 

Advective transport by g(t) and g{t) is defined by 

a(x^, t) = ao • g 1 ^) and a(x, i) — oq • <7 _1 (i) , 

where ao = a(xo,0) = a(xo,0), with a, a G V* and the 
factorization 5 (t) = E(t)-g(t) implies a(x, i) = a-S(x, t). 
Since a and a refer to the same initial conditions, ao, 



quantity. Here .F(x, t) is the tensor transformation 

factor of the advected quantity a under the change of 
variables 5 : x — > x'. For example, the density, D, 
transforms as 



D s det(TH)(x,£) = £>(x,£) 
and L> advects as 



.F(x, t) = det(TH) , (9) 
9 t £) = -div(Du) . (10) 



The transformation factors are 1, det(TS) and K = 
det(TS ) TS _1 , for an advected scalar, density and vec- 
tor field, respectively. In each case, the corresponding 
transformation factor T appears in a variational rela- 
tion for an advected quantity, expressed via equation (|^) 
as 



So, 6 = 6 (J- 1 ■ a) = T- 1 ■ 6a + {5!F~ 



(11) 



This formula will be instrumental in establishing the 
main result. 

Lagrangian averaged Euler-Poincare Theorem. 

Let the following list of assumptions hold ^|J|]. 

• There is a right representation of Lie group G on 
the vector space V and G acts in the natural way on the 
right on TG x V*: (v g ,a)h — (v g h,ah). 

• The function L : TG x V* -> R is right G- 
invariant. 

• In particular, if ao G V* , define the Lagrangian 
L aQ : TG —> R by L aa (v g ) = L(v g ,a ). Then L ao is right 
invariant under the lift to TG of the right action of G ao 
on G, where G a „ is the isotropy group of ao- 

• Right G-invariance of L permits one to define 
t : g x V* -> R by l{v g g~ x , OofT 1 ) = L(v g ,a ). Con- 
versely, this relation defines for any I : g x V* — > R a 
right G-invariant function L : TG x V" ^ R. 

• For a curve g(t) G G, let u(t) = g^g^)^ 1 G 
TG/G = g and define the curve a(t) as the unique so- 
lution of the linear differential equation with time de- 
pendent coefficients a(t) = —a(t)u(t) where the action 
of u G g on the initial condition a(0) — ao G V* is 
denoted by concatenation from the right. This solution 
can be written as the advective transport relation, 
a(t) = aog(t)~ l . 

• The GLM facto rization holds, g(t) = E(t) ■ g(t) 
with g(t) = S (t) • g(t) = g(t) and |(t) = 0(f). 

LAEP Theorem. 

The following are equivalent: 

i The averaged Hamilton's principle holds 



L ao (g(t),g(t))dt = 



(12) 



a 0' g — o( x i t) — a ■ H (x, i) = ^(x, i) -a^(x, i) . (8) satisfied on T* G ', 



for variations 5g{t) of g(t) vanishing at the endpoints. 
ii The averaged Euler-Lagrange equations for L ao a 



Note that the right side of this equation is potentially 
rapidly varying, but the left side is a mean advected 



SL a 



Sg 



■TE 



d 5L a 
dt Sq 



TS = 



(13) 
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iii The averaged constrained variational principle 
r*2 



S / l(u(t),a{t))dt = 
Jti 

holds, using variational relations of the form 

Su = (d t + ad u ) V ' + (TE • (d t + ad s )ry) S- 1 



(14) 



Sa 



ar l = 5(T~ 1 -a)- -1 



a r\ 



(15) 



where Lie derivatives of advected quantities by the vector 
fields r)'(t) = SE S , 77(7;) = Sgg^ 1 and 



V = $99 1 = V' + ( T - ■ V)- 1 , 



(16) 



are indicated by concatenation on the right and these 

three vector fields all vanish at the endpoints. 

iv The Euler-P oincare (EP) equation holds on g x V* 



< d ,*\ 81 st 
— + ad,. J — = — o a . 
Kdt V Su da 



(17) 



and the Lagrangian averaged Euler-P oincare (LAEP) 
equation holds on g x V* 



9 _L A* 

dt + a " 



o a . 



(18) 



Notation. In equations ( |l7j ) and (18), the operations 
ad* and o are defined by using the L^ pairing (/, g) = 
J fgd 3 x. The ad* operation is defined as (minus) the 
L2 dual of the Lie algebra operation, ad, or commu- 
tator, ad„r; = — [14,77], for vector fields, — (ad* /it, 77) = 
(fi, &d u r]). The diamond operation o is defined as (mi- 
nus) the L2 dual of the Lie derivative, namely, (boa , 77^ = 
— (6, £ r) a) = —(b,arj), where £ v a denotes the Lie 
derivative with respect to vector field 77 of the tensor a, 
and a and b are dual tensors. 

Proof of the LAEP Theorem. The equivalence of i 
and ii holds for any configuration manifold, so it holds 
again here. To compute the averaged Euler-Lagrange 
equation ( |l3|) , we use the following variational relation 

obtained from the composition of maps g(t) = 5 (t) -g(t), 
cf. the velocity relation (]l|), 

6g(t) = 6E (t) ■ g(t) + TE (t) ■ Sg(t) . (19) 

Hence, we find 



= 6 L ao {g(t),g(t))dt 
Ju 



' SL ao SL ao 

— . 6g + — . Sg )dt 



:(( 



5L ao d SL a 



/ 5g dt 

rt? / Tr 



•<5S(7j)J -gdt (20) 



SLa 



d SL a 
dt Sq 



TE -Sgdt. 



In the last equality, the first integral vanishes for any 
5 S, thus ensuring that the Euler-Lagrange equations are 
satisfied before averaging is applied. Vanishing of the 
second integral for arbitrary 5g then yields the averaged 
Euler-Lagrange equations (13). 

The equivalence of iii and iv in the LAEP theorem 
now follows by substituting the variations ( |l5|) into (H), 
and integrating by parts, to find 



r t2 r 2 1 si \ 1 SI \ 



. \ 51 SI \ , 

(St + ad J — o a , 77 ' ) dt 

Su 6a I 



TE 



(21) 



)),//. ,22) 



Thus, the independent variations 77' in (|T]) and 77 in 
result in the EP motion equation (|l7j ) and the LAEP 
motion equation (|Ts|), respectively. 

Finally we show that i and iii are equivalent. First 
note that the G-invariance of L : TG x V* — > R and the 
definition of a(t) = aoglt)^ 1 imply that the integrands in 
( |l2| ) and (|l4|) are equal. In fact, this holds both before 
and after averaging. Moreover, all variations 5g(t) G TG 
of g(t) with fixed endpoints induce and are induced by 
variations Su(t) S g of u(t) of the form Su = dn/dt+a,d u 77 
with r](t) £ g vanishing at the endpoints. The relation 
between Sg(t) and 77(7:) is given by rj(t) = Sg^g^)^ 1 . 
The corresponding statements also hold for the primc- 
and tilde-variables in the variational relations (|l5|) that 
are used in the calculation of the other equivalences. This 
finishes the proof of the LAEP Theorem. 
Lie derivative vs ad*. The equality ad*/i = £ u /i holds 
for any one-form density fi (such as /1 = St/ Su, the varia- 
tional derivative). Thus, the LAEP motion equation ( |lS| ) 
may be written equivalently using Lie derivatives as 



o a . 



In this notation, the advection of mass takes the form 
(dt + £u)D = and one immediately obtains the 
LA Kelvin-Noether Circulation Theorem. 



On defining v 



81 



TE)/D, 



d_ 

dt 



c(u) 



c(u) D \6at 



o a , 



for any closed curve c(u) following the LA fluid motion. 
The same result is obtained by applying LA directly to 
Kelvin's Theorem for the EP motion equation (Bl 
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Applying the LAEP theorem to incompressible 
fluids. The Lagrangian averaged Euler (LAE) equations 
for an incompressible fluid are derived from the LAEP 
theorem using the reduced averaged Lagrangian, 



t = 



d 3 3 



1 ^ 



D \u- + pt (d T'E D 



(23) 



The pressure constraint implies that the mean advected 
density is related to the mean fluid trajectory by D = 
detTS . Thus, in general, the LAE fluid velocity has a 
nonzero divergence [jfj, since ([l0|) for D implies 



div u 



T)\dt 



u L ■ V )D j£ 0. 



(24) 



In principle, one may restrict g(t) and both its factors 
H(i) and g(t) to the space of volume-preserving diffco- 
morphisms, for which detTS = 1. However, some LA 
processes may not respect this restriction and, in general, 
divu ^ 0. In lq j7|], though, divu = is accomplished. 
For £ in (|23j) the LAEP equation (fill) gives 



d _ 



d „ 



dxi % J dx % 



d 

dx l 



0. 



(25) 



with mean fluid quantities Vi and tt defined as 
1 Si FTZ^Ti . Si 1 

Jd 



When TS = Id 
finds v = u£ 



V£, for a vector field £ = — x, one 



-jjjrS.j'V^ = u — p . The term p is called 
the pseudomomentum (ij. See, e.g., for discussions 

of the role of pseudomomentum in GLM theory. 
Momentum balance. 

Following the EP theory of Holm, Marsden and Ratiu 
leads to momentum balance for the LAE equations, 



d 



d 



(bvi 



7,J 



■fS\ 



Dd\ut\ 
~2 dx i 



(26) 



exp 



where subscript exp refers to the explicit spatial depen- 
dence that yields a mean force arising from the S— terms 
in |u£| 2 = \D L 3 / Dt | 2 that appear in equation (Hj. 
Recent progress toward closure. Of course, the LAE 
equations ( p5| ) are not yet closed. As indicated in their 
momentum balance relation ([26]), they depend on the un- 
specified Lagrangian statistical properties appearing as 
the E!— terms in the definitions of v and tt. Until these 
properties are modeled or prescribed, the LAE equations 
are incomplete. 

Progress in formulating and analyzing a closed system 
of fluid equations related to the LAE equations has re- 
cently been made in the EP context. Such closed model 
LAE equations were first obtained in Holm, Marsden and 
Ratiu ^|J|] . For more discussion of this type of equation 
and its recent developments as a turbulence model, see 



papers by Chen et al. JlO) — |T~3|] , Shkoller JTJ] , Foias et al. 

Marsden, Ratiu and Shkoller f|] and Marsden 
and Shkoller 0. A self-consistent variant of the LAE 
closure was also introduced in Gjaja and Holm |^| in the 
Lagrangian fluid specification, see Holm [|]|]] for further 
discussion of that approach. 

The LAEP approach is versatile enough to derive LA 
equations for compressible fluid motion, as well. This 
was already shown in the original GLM theory For 
brevity, we only remark that the LAEP approach also 
preserves helicity conservation for compressible flows and 
it preserves magnetic helicity and cross-helicity conser- 
vation when applied to magnetohydrodynamics (MHD). 
For more details in this regard, see ||. 
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